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We study the domain of existence of a solution to a Riemann problem for the 
pressure gradient equation in two space dimensions. The Riemann problem is 
the expansion of a quadrant of gas of constant state into the other three vacuum 
quadrants. The global existence of a smooth solution was established in Dai and 
Zhang [Arch. Rational Mech. Anal., 155(2000), 277-298] up to the free boundary 
of vacuum. We prove that the vacuum boundary where the system is degenerate 
is the trivial coordinate axes. 
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1 Introduction 

The pressure gradient system 



>< 



{Ut+Px = 0, 
Vt+Py = 0, (1.1) 
E t + {pu) x + (pv) y = 0, 

where E = (u 2 + v 2 )/2 + p, appeared first in the flux-splitting method of Li and 
Cao [5] and Agarwal and Halt [1] in numerical computation of the Euler system of a 



* School of Mathematical Sciences, Fudan University, Shanghai 200433, P. R. China; School of 
Mathematics and Statistics, Northeast Normal University, Changchun 130024, P. R. China; Key 
Laboratory of Mathematics for Nonlinear Sciences (Fudan University), Ministry of Education, P. R. 
China, email: leizhn@yahoo.com. 

^Department of Mathematics, Pennsylvania State University, State college, PA 16802, email: 
yzheng@math.psu. edu. 



1 



compressible gas. Later, an asymptotic derivation was given in Zheng [TUl EE3] from the 
two-dimensional full Euler system for an ideal fluid 

lH + V-{pU) = 0. 

{pU) t + V-{pU®U + P I) = 0, 

( P E) t + V-(pEU + pU) = 0, 

where U = (u,v),E = {u 2 + v 2 )/2 + p/ ((7 — l)p), and 7 > 1 is a gas constant. We refer 
the reader to the books of Zheng [TT] and Li et. al. [7] for more background information, 
and the papers pUl El El [121 El EE] for recent studies. After being decoupled from the 
pressure gradient system (11.11) . the pressure satisfies the following second order quasi- 
linear hyperbolic equation 

(f) t - A c^ = °- ( L2 ) 

Dai and Zhang [2] studied a Riemann problem for system (11 .11) . see also Yang and 
Zhang |8] by the hodograph method. In the self similar variables £ = x/t,r] = y/t, the 
value of the pressure variable of the Riemann data is 

P{£,V) = £ 2 , fo r < f, 77 < Jpl, (f - VPi) 2 + V 2 =Pi, 

(1.3) 

P(£,V) = V 2 , for < f, rj < v /pT, £ 2 + (77 - v /pT) 2 = px. 

Here pi is any positive number. They showed that the Goursat problem for system 
(II .2p admits a global solution in the self-similar plane, which is smooth with a possible 
vacuum near the origin (see Figure 1, where a = tJpi)- We are interested in the size of 
the vacuum boundary {(£,77) | p(£,,v) = 0} where the pressure gradient system (11.21) is 
degenerate. Somewhat surprisingly, our result shows that the vacuum bubble is trivial 
and the entire vacuum boundary is the trivial coordinate axes in the self-similar plane 
(see Figure 2, where a = ^/Pl), which is stated in our main theorem at the end of 
Section 3. 

Further motivation for the study of the current problem is that the study of bound- 
aries such as a sonic curve is important in establishing the global existence of a solution 
to a general two dimensional Riemann problem of the pressure gradient system. In ad- 
dition, the solution of the current problem covers wave interaction problems in which 
only some fractions of the plane waves are involved. Wave interactions of these kinds 
are common in two dimensional Riemann problems. Finally, the study of the pressure 
gradient system has motivated work on two dimensional full Euler systems, see Li [1] 
and Zheng [MfT5] . 



2 Integration along characteristics 

In the self-similar variables £ = x/t, 77 = y/t, the pressure gradient equation ( 11.21) 
takes the form 

+ ydy) 2 P a m , (£dt + ydy)p ((^ + nd v )p) 2 

A (5,»?)P-^ 2 = U - V- 1 ) 

p p p z 
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Figure 1 . Solution with a presumed vacuum bubble. 

Figure 1: 
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Figure 2. The global solution. 

Figure 2: 
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In the polar coordinates 

9 = arctan |, 

equation (12.11) can be decomposed along the characteristics into the following form (see 

mm) 

d + d_p = mp r d^p, 

(2.2) 

d-d + p = —mp r d + p, 



provided that 
where 



p < r 2 , 



Ar 4 
m = — -. 

2p 2 ' 



and 




(2-3) 

r 2 (r 2 —p) ' 

Note that the equation is invariant under the following scaling transformation 
(£,»m0— 4=,-). (Pi>0) 

y/Pl y/Pl Pi 



Thus, without loss of generality, the corresponding boundary condition (11.31) in the 
polar coordinates can be set as the form 

p = £ 2 = r 2 cos 2 9 = 4 cos 4 9, on r = 2cos6», vr/4 < 9 < vr/2; 

(2.4) 

p = rf = r 2 sin 2 9 = 4 sin 4 9, on r = 2sin6>, < 9 < ix/A. 

The solution exists in the interaction zone up to a possible vacuum bubble. See Figure 
3. 

The characteristic form (12. 2p of the pressure gradient equation enjoys a number of 
useful properties. For example, the quantities d±p keep their positivities/negativities 
along characteristics of a plus/minus family, and the sign-persevering quantities yield 
monotonicity of the primary variable p (see [6], also [3] where the authors propose to 
call them Riemann sign-persevering variables), and the fact that a state adjacent to a 
constant state for the pressure gradient system must be a simple wave in which p is 
constant along the characteristics of a plus/minus family [B|. 

The characteristic decomposition (12.21) has played an important role and was a 
powerful tool for building the existence of smooth solutions in the work of Dai and 
Zhang [2j. We are interested in the size of the vacuum boundary {(r, 9) \ p(r,9) = 0} 
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Figure 3. Presumed vacuum in the polar coordinates. 

Figure 3: 
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where the pressure gradient equation is degenerate. For this purpose, we rewrite (j2.2p 

as 

d + d_p = qd + pd_p — g(<9_p) 2 , 
d-d + p = qd + pd-p — q(d + p) 2 , 



(2.5) 



where 



4p(r 2 — p) 

Define the characteristic curves rt{9) and r b + {9) as 



(2.6) 



drt (0) 
dd 



2 sin 9 a , 



and 



dr b + {9) 
d9 



r\(9 b ) 



i 

\(r b + (8),d) ' 

= 2 cos9 b . 



(2.7) 



We point out here that for convenience, we also use the notation r"(r, 0) (r^(r, #), 
respectively) which represents the characteristic passing through the point (r, 9) and 
intersecting the lower (upper, respectively) boundary at point a (b, respectively). See 
Figure 3. 

Now let us rewrite the system (I2.5P as the following form 



d +[w^ ex P Ie b 9<9+P( r +0)> <f)d(j) J = gexp qd + p(r b + ((j)), (p)dcj) 

= gexp Jg 6 * g<9__p(r" (0). 



Integrate the above equations along the positive and negative characteristics r b + {9) and 
r"(0) from 9 b and a to 9, respectively, with respect to 9, one can get the iterative 
expressions of d + p and d_p: 

( i 



d.p ex P So b <?<9+p(r+(0), = sb ( 2 cos ^' 
+ q {r b + , V>) exp g<9+p (r*. (0) , 0) d0dV 

ab ex P fe a qd - p ( r - (^) ' ^) rf< ^ = ab ( 2 sin e a,0 a ) 
. + Je%( r -(^)'^) ex P/et ^-P( r -(^)> 



(2- 



On the other hand, noting the boundary condition (j2.4j) . a straight forward calcu- 
lation shows that 



cxp 



qd + p(r b + ((p),4>)d(j) 



(2.9) 



exp^ 



1 

- + - 

p r 1 



V 



P*«(0),0) rl rM^ e ) 

V^COS^ I 4 J p (2cos9 b ,9 b ) 



d + p(r^(0),0)d0 

1 



r 2 (r^) — Pi 



dpi 



Similarly, one has 

exp / g<9_p(r^(0),0)d0 (2.10) 



p*(rt{0),9) rl fp{n(°)>°) 1 

ex P 1 7 / -?7-^ - d Pl 



V2sm6 a U J p{2 ^e a ,e a ) r 2 (r a _)-p 1 



Then, we have 



g(r^),V){exp qd + p{r h + {4>), 0)#}r# (2.11) 



r p 



and 



xexp \I /, . 

k 4 Jp(2cos6 b ,0 b ) r \ r +) Pi 



q(r a _(ip),?p)exp I g<9_p (r°_ (0) , 0) (20^ (2.12) 



1 / r p * 



4^2 sin a 7e a ^ 2 -P 



(r«(V>),V) 



x ex P 1 7 / -^r^ r#- 

U J p{2s in9 a ,e a ) r 2 (r a _)-p l J 

Finally, by substituting (I2.10p and (I2.12p into the second equality of (12. 8p . we arrive 
at a new iterative expression for d + p : 

d+p(r,9) (2.13) 

exp J e qd_p (r°_ (0) , 0) ri0 
^ (2 sin 6 a , tt ) + j [ q (rl ty) , if,) exp JjJ g«9_p (r« (0) , 0) ri0# 

4pi (r, 0) exp { J f p ^J a) ^jp^ dpi _ 



! sin tf a cos f 



+ £ (ri (VO , V>) exp { | f^r^) ^fcr^ } # 



Similarly, by substituting (12. 9p and (12. lip into the first equality of (12. 8p . one obtains 
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the new iterative formula for cLp : 



<9_p(r, 9) 

~ exp Jg b qd + p (r b + (0) , 0) d(j) 



(2.14) 




dpi \dip 



3 The boundary of the vacuum bubble 

In this section, we use the iterative formulas ( 12. 13f) and (12.141) to prove that the vacuum 
bubble {(r, 9) \ p(r, 9) = 0, 9 G (0, vr/2)} is in fact the trivial origin {(0, 0)} in the self- 
similar plane. We use the method of contradiction. Assume to the contrary that there 
is a bubble with boundary r (9) > for all 9 G (0, ~), and r (9) > for some 9 G (0, |). 
That means p(tq(9),9) = for 9 G (0, ~), and the solution is smooth in the domain 
bounded by the bubble r (9) and the upper and lower characteristic boundaries. We 
intend to deduce contradictions, which in turn proves that r (9) = for all 9 G (0, 7r/2) 
and the vacuum bubble {(tq(9), 9)} is in fact the trivial origin {(0, 0)} in the self-similar 
plane. 

Before we start the above procedure, we point out two observations which roughly 
imply the nonexistence of the vacuum bubble. For presenting the observations, we 
assume further that r (9) > for all 9 G (0, |). 

First, we can compute easily the characteristic slope 



along the upper characteristic boundary and consider 9 — > tt/2. Similar result holds 
on the lower boundary. Now assuming there is a bubble, noting that along the plus 
characteristics, by definition of A in (I2.3p . we see easily that 



as p — » (near the vacuum bubble) except for (0, 0) and (0, |). The above computation 
reveals that there is some kind of inconsistency for the slopes of the characteristics at 
(0, 0) and (0, ^) in the polar coordinate plane. 

Next, let us calculate the decay rate of p along the middle line 9 = tt/4. By using 
the symmetry of system (12. ip (see Figure 3), we can obtain that 




d±p ~ ±M p 1/2 
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asymptotically. To show the details, we propose 

d±p ~ ±M p^ +s , 

as (r, 9) tends to a point of the vacuum bubble. Then, by (12.131) . we have 

1 f i rp(r a -{6),e) i Wf i 

d+p(r, 9) = 4p4 ( r , 9) exp <^ - / — — dp x \ / <^ — 

•/ p(2 sin dafia) r v-)~Pi }l ^v^sin U a cosa a 

f e r 2 p-i rl fp( r -^^) 1 >, 

+ / U2_„\a - ( r -W^)exp - / . I2?^T— K 



9a (r 2 -p)d_p K ' Uy p(2sineaA) r 2 (r°) -p x 

M i+5 ,1-1 , , 

— -«4 T 4- high order terms, 
45 



which implies that 5 should be \. 
Thus, by using (12. 3p . we have 



asymptotically as (r, #) tends to a point of the vacuum bubble. Therefore, 

p ~ c exp I I (3.1) 

asymptotically as (r,9) tends to a point of the vacuum bubble on the line 9 = 7r/4, 
which implies that there is no interior vacuum at least at 9 — j. 

We point out incidentally that Zheng's previous numerical computation of the bub- 
ble, referred to in Dai and Zhang [2] , is probably caused by the fast exponential decay 

In what follows, we concentrate on establishing the above formal intuition rigor- 
ously. With the symmetry of system (12. 5p . let us restrict our arguments on 9 G (0, -]. 

Fix a point (r, 9) on the bubble boundary. Let us denote D(r, 9) the bounded 
domain determined by the positive and negative characteristic curves starting from 
(r,9) and (V2, f). See Figure 4. 

For < e < 1, define a curve r t (9), 9 G (0, |), by 

p(r 6 (0),0)=e. (3.2) 
From (I2.13P and (12. 14p . it is easy to see that 

'^ M)> °' for 0<9<*. (3.3) 
d^p(r, 9) < 0, ~ 4 

Thus, by (12.31) and (13.31) . we have 

Pr (r,9)>0, (3.4) 
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2 1/2 



Figure 4. Domain and notation 

Figure 4: 
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which implies that the curve r e (9), 9 e (0, |), defined in (13.21) is smooth if e G (0, 1). 
Here we point out that sometimes we still denote the characteristics passing through 
(r e (9),9) or any other point (r, 6) by r"(6 ) ) and r b + {6) or r"(r, 9) and r^_(r, 6 1 ) , when 
in fact a = a t is dependent on e, etc. We caution the reader that the intersection is 
determined by the point (r, 9) and the characteristic, which is hidden for notational 
convenience and always well understood. 

Now we fix eo = 1 (eo can be any fixed positive constant in (0, 1], which do not 
change any of the following arguments). Define 



where 



Then we have 



Mi = max {p 2 d + p, —p 2 <9 oj, 
over s 



S := {{r,9) | 9 E (0,7r/2), r £Q {9) <r<V2}. 



(3.5) 



d+p < Mip 
-d-p < M\p 



(3.6) 



for all (r, 9) with r > r eQ (9). Note that Mi depends only on e and does not depend on 
(7,9). 



Next, let 



A(r,9) 



B(r,9) 



cxp 



4 Jp(2aine a ,fl a ) r 2 (r a )_ pi 



dpi | 



ex P 1 4 J P (2 cos e^e,,) r 2 (r (> )_ Pl api 



(3.7) 



2^ 



cos^ f;, sin tf(, 



Then, let 



. 4p-i(2sin0 o A) 4p-3(2cos^,^) 1 
M 2 = max -j max _ jt — — , max _ — — — 1 , 



A{r,0) 



(r,6)eB(r,0) 



Br, 



and 



M 3 = max {Mi,M 2 + l}. 



(3. 



(3.9) 



Just as what we had pointed out before, the intersections a and b in (13.71) and (13. 8p 
vary with (r, 9) and characteristics, which are not expressed explicitly for notational 
convenience. We intend to prove that the inequalities (13.61) are still valid for all points 
(r, 9) G D(r, 9) with M\ being replaced by M3. Namely, there hold 



d + p < M 3 p2. 
—d_p < M 3 p2. 



(3.10) 



Note that the positive constant M 3 (> Mi) is independent of (r,9) 6 D(r, 0), but 
depends on the fixed (r, #). 
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We start to prove (I3.1UI) . Suppose that (I3.10p is correct up to a line segment 
r e (6) C D(r, 6), then we improve f)3.10p to strict inequalities on the line segment r e (8) 
in the same domain. In fact, for (r, 9) e D(r, 6) with r < r eo (8), let us compute 

<9 +P (r,0) = 4p*(r,0)/{A(r,0)+ / ^-i(r°(V0, V>) (3.11) 
= 4 P i (r, 6)/{A(r, 6) + f ^I'^X-p) (r ° W,) ' ^ 



where 



/(0" 



x ex P <! 7 / , >. -^77 d Pi 



4 



Er 2 



< M 3 p*(r,e)/{ 4 l ; + / X^^W^) 



p(2sin0 a ,0 a ) r P 
X eX P l 7 / / X 2f rf Pl \ d P 

M 3 pi(r,0) 

a^M) _i -I J 

4 +il c/ J Jp(2sme a ,e a ) aP 4a P 

M 3 pi(r,0) 



M 3 A(r,6>) 
4 



+ /(0-)[p-3(r,0)-p-i(2sin0 o ,0 a 



r 2 fl /■ p ( r - (e_),e ~) 



exp{ i f_ ^)^*.} (3-12) 



r 2 - p v ' 1 l 4 J p (r«_(e),e) 

2 r P (r a _(e-),e-) 1 

> ex P i 7 / , v ^2 rf Pl 

/p(r»(0),0j V2 
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-p(r«(0-),0-)-p(r^(0),0) 
exp < — ^ > 1 



with some 6 a < 6~ < 6 and A(r, 0) is defined in (1377j) . Thus, by Q), fIBTTTD and (15321) . 
we have 



M 3P 4( r ,0) 
+ [p-*(r,0)-p-*(2sin0 a ,0 a )] 

i 



w) < M ,,„, 8 , . , ~rr"'., (3.13) 
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Similarly, by O and (l2~TTj) . we have 

i>e 2 -2 

-d-p{r,9) = 4p 1 *(r,6)/\B(r,9)- L^-L^^)^) (3.14) 



< 



— p 

dpi\di/j\ 

)-Pi > > 

M 3 pi(r,9) 



M 3 B(r,B) 
4 



where 



+ g(0+)[p-?(r,e)-p-*(2coB9 b ,e l 

r 2 rl r p ( r +( e+) > e+ ) 1 

fl r P (r b + (e+),e+) 1 

> ex P 1 7 / , v ^2 rf Pl 

f pK(fl+),g + )-p(r^),fl) 
= exp < — ^- \ > 1 



with some 9 < 9 + < 9^ and B(r,9) is given in (13. 7p . Thus, the similar estimate as 
(13TT3]) holds: 

-0_p(r,0) < M 3 p^(r,9). (3.15) 

Since M 3 depends only on (f,9), particularly, is independent of (r,9) G B(r, #), the 
proof of (13TTUD follows by fl37T3|) and fl3~T5|) . 

Now with the aid of (12.31) . we add up (13. 13j) and (13.151) to yield 

M 3 2M 3 

Pr < , =P < —Z2~P, 

yr 2 (r 2 — p) t 



for all (r,9) in a small neighborhood of (r,9) in B)(r,9). Thus, a simple integration of 
the above inequality with respect to r from r to r yields 

M 3 

p(r, 9) > p(r, 9) exp{ — -^(r — r)} for r < r. (3.16) 

On the other hand, by ( 13.41) and the fact that p(ro(9), 9) = 0, we have 

p(r, 0) > 

for all (r, 0) with r > r (6'), which together with (13.161) result in that p(f, 9) > 0. Since 
(r, 9) is a point on the bubble, we arrive at a contradiction. 
Summing up, we in fact proved the following theorem. 

Theorem 3.1. The Riemann problem ( 11.2ft . ( 12.4ft /or t/ie pressure gradient equation 
admits a unique smooth solution. The pressure of the solution is strictly positive in 
£ > 0,r? > 0. 
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